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Figure 1: The FOSH operator shown acting on the connection coefficients. Dot indicates zero. 
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Abstract 



By adopting Nester's 4— dimensional special orthonormal frames, the tetrad equa- 
tions for vacuum gravity are put into explicitly causal and symmetric hyperbolic 
form, independent of any time slicing or other gauge or coordinate specialization. 



0.1. Introduction 



We have previously given a well set and causal exterior differential system for 
vacuum gravity, generated by a closed set of differential forms describing the 
immersion of 4— dimensional spacetime into a flat 10 dimensional space PQ. The 
orthonormal frame bundle of the latter has a canonical basis of 10 (translation) 
1— forms uj^ and 45 (rotation) 1— forms u^, satisfying the structure equations 
of the Lie group ISO( 10). Dividing the range fx, v = 1, • • • , 10 into two ranges 
i, j = 1, • • • , 4 and A, B = 5, ■ • • , 10, these are 



du l + u) 


Au j 


+ u B Au B 


= 


du A + uf 


Aw 3 


+ u B A uj b 


= 


du\ + (J* 


a 4 


+ u B A u B 


= 


du l A + Uj 
due + uf 


A u{ 

A 4 


+ uj b A uj b 
+ u B A uj c 


= 
= 



The exterior differential system is generated by 6 (immersion) 1— forms uj a , 
their closure 2— forms du A = —uf A uj\ and 4 closed 3— forms ensuring Ricci- 
flatness, namely, Rij A u>k£ ljkl , where = dtUj + u\ A u k , — —u l A A uof defines 
the Riemann 2— forms. 

The significant result is the calculation of the Cartan characteristic integers 

s = {s , si, s 2 , s 3 } = {6, 6, 10, 8}. 

This shows the solutions to be 25 dimensional, regular (i.e., in principle con- 
structed from a series of Cauchy-Kowalewsky integrations) and causal (i.e., = 0, 
so the solutions are determined from suitable data set on 3 dimensional immersed 
manifolds). The solutions are involutory with respect to u l ,Uj and u B : these 
forms remain independent when pulled back into a solution manifold and can be 
adopted as a basis there. 

The six basis forms uj % - and 15 basis forms u B , which occur of course in the 
structure relations, do not appear explicitly in the exterior differential system, 
showing the solutions to be bundles having 21 dimensional fibers over a four 
dimensional base. Evidently this expresses arbitrary 0(4) (for Lorentzian, 0(3, 1)) 
rotations of the tetrad frame u/ and 0(6) rotations of the immersion co- frame u A 
at each point of the base. On a four dimensional cross-section all the forms can 
be expanded on the of basis, the u l j being a metric connection. 
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In Section 2 we give a new immersion exterior differential system for vac- 
uum gravity that incorporates the higher dimensional special orthonormal frame 
(HSOF) conditions proposed by Nester [2] as a generalization of the special or- 
thonormal frame conditions (SOF) in three dimensional Riemannian geometry [3| 
[I] Calculation of the Cartan characteristic integers for this exte- 
rior differential system for 4 dimensional Riemannian geometry shows 
that it is also well set and causal, so that, as Nester conjectured, such 
frames can be imposed without impediment in extended regions of vac- 
uum spacetime. Solutions are now 19 dimensional, the fibers expressing only 
arbitrary 0(6) rotations of the co-frames. 

These exterior differential systems include all integrability conditions for the 
determination of a metric on the base, or on any 4 dimensional cross section. 
In terms of base space coordinates x^ (fi,f = 1,2,3,4), invertible matrices of 
functions l X ft (x) and iX^(x) exist such that *A MJ A^ = 5*-, = 1,2,3,4). We 
introduce the Minkowski metric, i.e., rfi = diag(l, 1, 1, — 1), to raise and lower 
left, or Lorentz, indices. Then the metric is given by g^ u = r) kl ^X^ iX u . Inserting 
uf = dx^ in the structure equations and in the exterior differential system 
gives the partial differential equations for coordinate components of the tetrad 
field. 

In Section 3 we reformulate the higher dimensional special orthonormal frame 
system in orthonormal tetrad components, using rfK We use the dyadic formalism 
of references jH] and [7j, in which the unit 4— vector field A X^ is given a special 
meaning: it traces a congruence of timelike world lines, a 3 parameter "fluid" of 
point observers, to which physical interpretations of the 24 components of the 
connection are attributed. The three spacelike unit vectors a X fl (a = 1,2,3) at 
each point complete a local orthonormal frame for the observer there. We expand 
the six 1— forms u\ on the u a , uj a basis (a, b = 1,2,3), and we expand the six 
2— forms Rj on the u a A u> b , u a A a; 4 basis, subject to the conditions for Ricci- 
flatness. 

The 24 connection components are grouped into the following ensembles: 3x3 
dyadics K and N, and 3— vectors a and u. The dyadic K has components K ab , and 
can be resolved into 

where S a b is the symmetric rate-of-strain 3— tensor of the observer fluid, and fl c 
is its axial vector of vorticity. Or, in dyadic notation 

k = s-nxi, 
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where I is the unit dyadic. The dyadic N is formed from the nine spacelike Ricci 
rotation coefficients of the ui a basis. It has components N ao and can similarly be 
resolved into a symmetric part and an axial vector n c 

Nab = N s a l m - n c e acb . 

Again we write in dyadic notation, 

N = N ssm -nx I. 

The vector uj has components u a and is the time-dependent angular velocity of the 
triad seen by an observer moving along uj 4 , with respect to a Fermi-propagated 
frame. Since u is a standard notation for angular velocity this usage should in 
context not be difficult to distinguish from the 1— forms ui a we have used up to 
this point. (When necessary, we can write the components of the angular velocity 
as uj a .) The vector a has components a a and is the acceleration of the point 
observers, i.e., their departure from geodesic motion. The vectors a and u are in 
principle determined operationally by the point observers, using spring balances 
and supported spinning particles in the local frames. The ten components of 
the Weyl tensor yield symmetric tracefree dyadics A and B, the "electric" and 
"magnetic" tidal fields. The quantities A, B, K, N, a, and uj are defined in terms 
of the iX^ u in j^j. 

The dyadic formalism is completed by the use of inner and outer 3— dimensional 
multiplication ( • and x ) , and by convective derivatives (also known as unit deriva- 
tions) in the timelike and spacelike directions ( 'and D ). Use of spatial covariant 
differentiation V, related to D by the spatial connection N, often is more efficient 
than the use of D. Detailed expositions of this formalism, including the various 
3— vector and dyadic relations, the relations of D and V (involving N), and the 
52 general first order dyadic differential equations for vacuum 4— geometry, are to 
be found in references [H] and j7]. 

There are some notational changes of which the reader should be aware. The 
dyadic we now (and in [J]) denote by N was in denoted by N*. For the dyadic 
N in [6 (which was symmetric) one should now understand 

N — (TV N) I + n x I 

and for the vector L in |J| one should now write n. 

Section 3 also gives the 12 new equations specializing to 4— dimensional special 
orthonormal frames. Six of the additional equations involve a and uj, which had 
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not otherwise appeared in the tetrad equations. This is also the case with the 
six conditions for the "Lorentzian frame" specialization (based on time slicing) 
discussed by van Putten and Eardley jH]. When time evolution of all dependent 
variables is explicit the causal structure of the system can be seen as resulting from 
local wave propagation subject to constraints expressing integrability conditions. 
The integrability properties of these sets of dyadic equations are briefly discussed 
in Section 3 in terms of Cartan's reduced characters and Cartan's test. 

In Section 4 we present the complete set of 64 dyadic equations resulting from 
substitution of dyadic components in the structure equations and the HSOF exte- 
rior differential system of Section 2. There are 34 equations involving the timelike 
derivatives of the 24 connection coefficients and the 10 Weyl components, together 
with an additional 30 transverse equations in which no time derivatives appear. 
By forming appropriate linear combinations, the final equations have been ar- 
ranged to show that they have first order symmetric hyperbolic (FOSH) structure 
similar to that of recent formalisms involving preferred time slices, designed for 
application to numerical gravity [HI] [H] [121 P^l - 

A number of known and new "hyperbolic reductions" of the vacuum field 
equations, and the choices of gauge they allow, have been surveyed by Friedrich 
[T4*j . He uses both tetrad formalism, and ADM variables based on preferred time 
slicing. The role of the Bianchi identities is emphasized: the equations for A (in 
[T4] . E) and B, which propagate as spin— 2 massless fields (cf., e.g., [Hj), are given 
in FOSH form. The FOSH equations obtained in the present paper necessarily 
include these. Nester's conditions on the tetrad components of the connection 
however leave no further freedom in choice of gauge (or frame), and seem not to 
have been previously used. They result in FOSH structure with constant 
coefficients, and force all the dependent variables to evolve along null 
cones. 

An analysis using Cartan's test shows that in the present case the observer 
fluid necessarily has vorticity, i.e., f2 ^ 0, so it is not 3— space orthogonal, and 
therefore Nester's special orthonormal frames cannot be based on a preferred slic- 
ing. Nevertheless, the relative simplicity of these equations may be advantageous. 
No lapse or shift variables or higher order derivatives have been introduced. Of 
course, useful related coordinates can still be found. Two are suggested imme- 
diately by the closed 2— forms in the exterior differential system, and we note 
in Appendix 1 that the dyadic conditions for a harmonic timelike coordinate are 
again of first order symmetric hyperbolic (FOSH) form, so can simply be added 
to our results. Further adopting three spacelike coordinates co-moving with the 
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point observers then allows an explicit line element to be written, given a solution 
of the FOSH equations. 

A similar application of special orthonormal frames can be made for 2 + 1 
gravity, again leading to a constant coefficient FOSH system. This is outlined in 
Appendix 2. 

0.2. Higher dimensional special orthonormal frames 

Using eight new variables jji and {% = 1, 2, 3,4), we prolong the previously given 
immersion exterior differential system for vacuum gravity, with two additional 
closed 2— forms, two additional 3— forms, and their closure 4— forms. The 3— forms 
and 4— forms essentially define the and Z{ in terms of the connection forms 
(the Uj pulled back into the solutions), and the 2— forms require them to satisfy 
Dirac type partial differential equations. When expanded in tetrad components 
in Section 3 it can be verified that this is precisely Nester's prescription. The 
generators of this exterior differential system in 63 dimensions are 

uf A J 
(dyi - uJiVj) A u* 
(dzi — ujj Zj ) A UJ 1 
(u4 A u k A J + \vi J A uo k A u l )e* jkl 
Uij Aw'A LU j - \ziUJ j A uj k A u) l e\ jkl 

(-u j s Aul Au k Au l + 2u{ Au k s Au s Auj 1 + ldyiAuo j Auo k AJ -2yiulAu s Auo k AJ)ei jkl 
-UisAujj Auj 1 Auj 3 +2uJij Auj 1 s Auj s Auj 3 -(\dziAuj j Au k Au l -ZiU J s Au s Au; k Au l )e l . jkl . (1) 
Monte Carlo calculations |T] of the Cartan characteristic integers yield 

8 = {6, 8, 14, 16}, 

therefore the exterior differential system is well set and causal; solutions are 19 
dimensional, fibered over 4 dimensions. The do not appear in the exterior 
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differential system, so the fibers express co-frame rotation. But the uA are explic- 
itly present so the frames are specialized and determined up to a simultaneous 
rotation at every point. 

If the two 2— forms are exact, two further variables can be added, say C and 
rj, together with 1— forms 

d( - y { uj\ 
drj — ZiUj 1 . 

These should be useful for introducing intrinsic coordinates into the HSOF for- 
mulation. 

0.3. The dyadic components of the connection forms Uj and 
Riemann forms R l - 

We expand the Uj on u\ defining the 24 3— dyadic and 3— vector components of K, 
N, u, and a, summarized in the Introduction and described in detail in references 
© and [Zj 

u) 12 = -U21 = N^u 1 + N 23 u 2 + N 33 u 3 + ufu 4 

<^23 = -^32 = N 2 lU) 2 + N 31 UJ 3 + NiiUJ 1 + Uj 1 ^ 4 

uj 31 = -u 13 = N 32 lj 3 + N12UJ 1 + N22U 2 + uj 2 uj 4 

Uu = — W41 = —Knu 1 — K 21 co 2 — K 3 iu 3 + a l uo 4 

u 2 4 = -w 42 = -K22L0 2 - K 32 to 3 - Ki 2 uJ l + a 2 u 4 

W34 = -W43 = -K 33 u 3 - K^u 1 - K 23 uj 2 + a 3 u 4 . (2) 

We also expand the Riemann 2— forms on a basis consisting of u> a A u 4 and u a A 
to (a, b = 1, 2, 3), to define (in the Ricci-flat case) Weyl dyadics A and B 

\R 2 i = B 31 u x Au 4 + B 32 co 2 Au 4 + B 33 u 3 Au 4 + A 31 u 2 Au 3 + A 32 uj 3 Au 1 + A 33 u x Aco 2 

\R 32 = B 12 u 2 Auj 4 + B 13 uj 3 Au 4 + B^u 1 Au 4 + A 12 u 3 Au l + A 13 u l Au 2 + A xl u 2 Au 3 
\R X3 = B 23 u 3 Auj 4 + B21UJ 1 Au 4 + B 22 u 2 Au 4 + A 23 u l Au 2 + A21U 2 Auj 3 + A 22 0J 3 Auo 1 
\R A1 = -A^u 1 Au 4 -A 12 u 2 Au 4 -A 13 u 3 Au 4 + B xl u 2 Au 3 + B 12 u 3 A^ + B^u 1 Au 2 
\Ra2 = -A22U 2 Auj 4 -A 23 u 3 Auj 4 -A2iu 4 Au 4 + B 22 uj 3 Au 1 + B 23 lo 1 Auj 2 + B 2 iuj 2 Alu 3 
I# 43 = -A^A^-A^A^-A^A^+B^A^+B^A^+B^Au 1 . (3) 
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The Riemann 2— forms given above are antisymmetric under interchange of their 
indices. The 10 Weyl components satisfy A a b = Af, a , A aa = 0, B a b = B^a, and 
B aa = 0. It can be verified that the Riemann symmetries are satisfied, namely, 

and 

R 4 b Auj b = 0, 

and also that the four 3— form conditions for Ricci-flatness, namely, 

R)Au k e j m = 

have been imposed. 

From the 3— forms of (1) we can now find the component expansions of the 
fields yi and z,, on an integral manifold of the exterior differential system. These 
are given by the following 

^ = (y a ,-y 4 ) = (Nxl-a,+TrK) (4) 

Zi = (z a ,-z 4 ) = (-KxI + cj,7VN). (5) 

We could alternatively have used vectors 2n = Nxl and 2f2 = Kxl. The twelve 
first order dyadic equations arising from the two additional HSOF 2— forms then 
are 

(Nxl-a)"+ V(TrK) + (TrK)a+ K-(Nxl-a) + wx (Nxl-a) =0 (6) 

V x (Nxl— a) + (TV K) Kxl = (7) 

(Kxl - uj)'- V (TV N) - (TrN) a + K- (Kxl -co) + uj x (Kxl) =0 (8) 

V x (Kxl -u) - (TrN) Kxl =0. (9) 

The sets of first order equations we have derived, where the dyadic components 
are taken as dependent variables and the four uo l form a basis, may alternatively 
be analyzed by another of Cartan's techniques. We must know that the set 
includes all integrability conditions, and it must also be possible to write the 
exterior differential system such that the left hand sides linearly involve exterior 
derivatives of the dependent variables (e.g., no terms of the form dK ab A dN cd are 
allowed), while the right hand sides only involve forms in the adopted basis (here 
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uj 1 , uj 1 Auji, etc.). So-called reduced characters s' { are then conveniently computed 
from the left-hand sides alone, and Cartan's test is to calculate 

If h is equal to the number of independent first order equations, this establishes 
involutivity, i.e., the well set nature of the problem. Moreover, if 

3-1 

s' g = n - g - ]T s[ = 0, 

i=0 

the Cauchy-Kowalewsky solutions are causal, i.e., determined by data given on a 
g — 1 dimensional surface. In the present context g = 4. 

The 52 general vacuum dyadic equations were of this form [S| [T^]. They result 
from computing the exterior derivatives of (2) and (3), so have six 2— forms and 
six 3— forms in 34 dependent variables. Their left hand sides are 

dN ba A u h - dud a A lu 4 

dK ba A u b + da a A 

dA ab A u b A uj a - | dB ab A u c A u d e h cd 

dB ab Au b A + | cL4 ab A a; c A u d e b cd . 

We calculate s' = {0, 6, 12, 10} , so h = 52. There are however s'^ = 6 arbitrary 
functions in the solution, and this system is therefore not causal. The new HSOF 
system of dyadic equations adjoins 12 equations in two additional 2— forms to 
express equations (6)-(9). Their left hand sides are 

d(2n a - a a ) Au a + d(Tr K) A c/ 

d(2Vl a - u a ) Au a - d(Tr N) A uj 4 . 

The reduced characters are s' = {0,8, 14, 12}, h = 64, s'^ = 0, and therefore the 
final system is causal. 
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0.4. The FOSH dyadic equations 



Linear combinations of the 52 general dyadic equations and the 12 new conditions 
due to Nester can now easily be made to put the result in FOSH form. The result 
is 34 equations involving the time derivatives and symmetric space derivatives of 
the dyadic components and A and B, and 30 "constraint" or transverse relations 
not involving time derivatives. We give them in the following, written in full with 
their quadratic right hand sides. For the connection components we obtain 

a-V • K T +V xw = -KxN-w-N + (Tr N)u; + 2K-n (10) 

Cj+V ■ N T -Vx a = -(Tr N)a+ (2Sl-u) • K-2(Tr K)fi - 2n-N (11) 

K - aV-2Vn + 2nV = -K • K — 2(Tr K)f2xl — uj xK + Kxw + aa-A (12) 

N + wV + 2 VfJ - 2fJ V = -K-N-u;xN + Kxa - auj - 2(Tr N)fix I + B. (13) 

We mention in passing that using D instead of V usually makes the right hand 
sides of these equations less concise. The prominent exception to this statement 
is the first equation, which becomes homogeneous and linear, namely, 

a-D-K r + Dxw=0. 

It has been convenient to use n and f2, and V can operate from the right as 
well as from the left, to express transposed indices. It is best to return to the 
nine components of K and of N to see the FOSH structure; e.g. write (2Vn) n = 
ViAT 23 - ViAT 32 , (2nV) 23 = V 3 ^3i - V 3 AT 13 , (2£W) 22 = V 2 K 31 - V 2 K 13 , etc. 
The left hand sides of the FOSH equations (10)-(13) are shown in Figure 1 as the 
product of a linear operator and a column vector. 
The 24 constraint equations are 

VxK=2fia-B (14) 

V x N = -|N T X N-|K X K + (fi-K) x \-2flu- A (15) 

V x (a-2n) =2(Tr K)fi (16) 

V x (w-2n) = -2(Tr N)0 (17) 

By taking a second time derivative of the FOSH equations, interchanging space 
and time derivatives, and substituting back both the FOSH and the constraint 
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equations, it can be seen that the dyadic variables all propagate causally along 
null cones. 

The 10 FOSH equations for the Weyl components are those for traceless trans- 
verse massless spin— 2 fields (dyadic and Bianchi identities, linearly combined) and 
are as follows 

2B-VxA + AxV = 
2Bxu-2ux B-2Axa + 2axA + K T -B + B-K-2(TrK)B + K X B + B X K 

X X 

(18) 

2A + V x B-B x V = 
2A x uj-2uj x A + 2Bxa-2axB + K T -A + A- K-2(TrK)A+ K X A + A X K 

X X 

(19) 

The left hand sides of the Bianchi equations are shown in Figure 2 as the product 
of a linear operator and a column vector. 
The final 6 constraint equations are 

V • A = -KxB -4fi-B (20) 

V • B = KxA + 4ft-A. (21) 

Finally, it should be remarked that our derivations from well set exterior differ- 
ential systems obviate any need to verify that the transverse (or constraint) equa- 
tions, i.e., those not involving time derivatives, are compatible with the FOSH 
system, and that they are propagated invariantly by it. 

0.5. Appendix 1: Harmonic and co-moving coordinates and 
an explicit line element 

The formulation of this paper is coordinate free and gauge independent. It may 
however also be of use to briefly record how a harmonic time coordinate, that is, 
one which satisfies a wave equation, and co-moving (with 4 A^) spacelike coordi- 
nates can be adopted. 

In coordinate language put t-^ v g^ u — (fi, v — 1, • • • , 4). Introducing fields 
(f) — t, A — Vt, the first order tetrad equations for this are 

A-V0=-K-A-wxA + 0a 
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4> - V • A = a • A - (TV K) 
V x A = 2<pQ. 

This set consists of four FOSH equations plus three constraints. It can be added 
to, and solved simultaneously with, the equations in Section 4. 

Co-moving spacelike coordinates, say x a (a — 1, • • • , 3), such that x a = are 
found by setting 

e a = \7x a , 

and the integrability conditions for this are 

e a = -K • e a - to x e a 



Vxe" = 0. 

The e a are, understandably, the first dependent variables we have found whose 
causal propagation is strictly timelike and not along the local null cone. We 
introduce coordinate components A-e a = A a and h al3 = e a -e 13 , and also calculate 
the inverse h a p {h ai h af3 = 5J). Then <fi, A a = h a pA^ and h a p, functions of x a , 
and t, enter the final line element: 

ds 2 = -4>- 2 dt 2 + 2<p- 2 A a dx a dt + (h aP - 4>~ 2 A a A p ) dx a dx p . 

In covariant 4— vector terms, the coordinate components of the tetrad vectors take 
the form 

' a\ a aA \ 
6 



a K o 



The orthonormal triad components of the vector A are written a X a A a . 



0.6. Appendix 2: 2 + 1 dimensional gravity 

An entirely parallel development can be made using special orthonormal frames 
(SOF) in 2+1 gravity. The exterior differential system for immersion of 3— dimensional 
flat spaces in the 21 dimensional orthonormal frame bundle over 6 dimensional flat 
space, which is generated by the usual immersion forms uj a and duj A {A = 4, 5, 6), 
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and the Riemann 2— forms R b (a,b,c = 1,2,3), has s = {3,6,3} and g = 9. We 
add a 2— form 

u b a A u c e a bc + y a uj b A co c e a bc , 

a 3— form 

u ab A/Aw'- zu 1 Aid 2 A u 3 , 

and their closures, to introduce 4 new variables y a and z, so the exterior differential 
system is prolonged to a total of 25 dimensions. The SOF conditions are imposed 
as two additional closed exterior forms, namely, 

dz 

d{y a u a ). 

These adjoin solutions of a linear Dirac equation j4J. Now the Cartan characters 
show the exterior differential system to be well set and causal with s = {4, 8, 7}, 
and g = 6. Solutions are 3— parameter co-frame rotation bundles over 3— space, 
determined by causal integration from 2— spaces. 

Orthonormal triad formalism for 3— dimensional geometry yields 9 equations 
for 9 Ricci rotation coefficients grouped in a dyadic N: 

V x N + ±N T X N =0. 

2 x 

The SOF conditions require the identification of z as Tr N and the y a as the 
components of n = |Nxl. To the nine equations for a flat 3— space are added six 
SOF equations: 

V(TrN) =0 
Vx n = 0. 

Inserting a —1 corresponding to timelike nature of the 3 direction, the complete 
set falls into FOSH form with 6 transverse equations. The left hand sides of the 
FOSH equations are shown in Figure 3 as the product of a linear operator and a 
column vector. 
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